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ABSTRACT 

The underlying mechanisms driving the quenching of dwarf-mass satellite galaxies 
remain poorly constrained, but recent studies suggest they are particularly inefficient 
for those satellites with stellar mass 10^ Mq. We investigate the characteristic evolu¬ 
tion of these systems with chemodynamical simulations and idealised models of their 
tidal/hydrodynamic interactions within the Mq group-mass hosts in which 

they are preferentially quenched. Our fiducial simulations highlight the role played by 
secular star formation and stellar bars, and demonstrate a transition from a gas-rich 
to passive, H I-deficient state (i.e. ASFR^-1, defHi^O.5) within 6 Gyr of first infall. 
Furthermore, in the 8-10 Gyr in which these systems have typically been resident 
within group hosts, the bulge-to-total ratio of an initially bulgeless disc can increase 
to 0.3<B/T<0.4, its specific angular momentum Ar reduce to ~0.5, and strong bisym¬ 
metries formed. Ultimately, this scenario yields satellites resembling dwarf SOs, a result 
that holds for a variety of infall inclinations/harassments albeit with broad scatter. The 
key assumptions here lie in the rapid removal of the satellite’s gaseous halo upon virial 
infall, and the satellite’s local intra-group medium density being defined by the host’s 
spherically-averaged profile. We demonstrate how quenching can be greatly enhanced 
if the satellite lies in an overdensity, consistent with recent cosmological-scale simula¬ 
tions but contrasting with observationally-inferred quenching mechanisms/timescales; 
an appraisal of these results with respect to the apparent preferential formation of 
dSOs/SOs in groups is also given. 
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1 INTRODUCTION 

Besides their initial and often rapid formation, the evo¬ 
lution of most galaxies can be defined by a gradual 
suppression of star formation (quenching), usually co¬ 
incident with a transformation from a disc to ellipti¬ 
cal morphology, much of which is observed to occur 
within within cosmologically overdense regions like groups 
and clusters (Butcher & Oemler 1984; Dressier et al. 1997; 
Lewis et al. 2002; Kauffmann et al. 2003a; Balogh et al. 
2004; Weinmann et al. 2006). 

In the context of the ACDM framework, these obser¬ 
vations partially reflect the early hierarchical assembly of 
massive galaxies via merging and their quenching via radio¬ 
mode AGN feedback (Balogh et al. 2004; Bluck et al. 2014). 
This mass-quenching appears to dominate at high-^ and for 
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masses above a characteristic 10^° M© (Peng et al. 2012; 
Wetzel et al. 2013). The (dwarf) galaxies that lie below this 
threshold are rarely found quenched in the field (Geha et al. 
2012), but tracing their recent quenching/transformation to 
specific mechanisms as satellites remains difficult. 

The main mechanisms are speculated to include tidal 
shock/heating by a host potential well (Mayer et al. 2001; 
Kazantzidis et al. 2011), harassment by successive tidal 
encounters (Moore, Lake & Katz 1998; Wetzel & White 
2010; Bekki & Couch 2011), and hydrodynamical inter¬ 
actions with the hot gas halo of a host (Intra-cluster 
medium; ICM). The latter, often invoked to explain 
a suppression of star formation (SF), include in par¬ 
ticular the ram pressure stripping of a satellites’ cold 
gas disc (Gunn & Gott 1972; Koopmann & Kenney 2004; 
Chung et al. 2007) and the stripping of the satellites’ hot 
gaseous corona that would otherwise fuel star formation 
(strangulation; Larson, Tinsley & Caldwell 1980). 

The efficiency of these processes within a given environ- 
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merit is expected to vary with properties of the constituent 
baryons (i.e. satellite velocity dispersion, ICM tempera¬ 
ture) and host halo. In practice, a bimodal (red and blue) 
galaxy population is observed in most overdense regions 
(van den Bosch et ah 2008; Cibinel et ah 2013; Bluck et ah 

2014) since z = 0.2 — 0.5 (McGee et al. 2011; Popesso et al. 
2012). This bimodality extends down to the host mass scale 
of 10^^ M©, in which the majority of all galaxies reside 
(Treu et al. 2000; Eke et al. 2004; Rasmussen et al. 2012). 

Recently, statistical analyses of large cosmological simu¬ 
lations and surveys (e.g. SDSS) have been adopted to estab¬ 
lish i) the timescale of satellite quenching/transformation, 
with respect to ii) the time since accretion to an overdense 
region (e.g. De Lucia et al. 2012). The finding that these 
timescales lie on the order of several Gyr is consistent with a 
slow quenching facilitated by strangulation (McGarthy et al. 
2008; Weinmann et al. 2010), as contrasted with the rela¬ 
tively rapid action of ram pressure stripping, which in some 
models produces too many passive satellites (Wang et al. 
2007). 

Other studies have gone further to suggest that the 
quenching timescale is independent of host halo masses 
above >10^^“^^ M© (e.g. Wetzel et al. 2013), where any 
apparent variability corresponds instead to the hierarchi¬ 
cal growth of such environments. In the event of this 
pre-processing, previously quenched satellites may be ac¬ 
creted as part of their host to a more massive host 
(Zabludoff & Mulchaey 1998; McGee et al. 2009; Peng et al. 
2012). This remains far from conclusive, however, with nu¬ 
merical studies differing in their prediction for the onset of 
these quenching mechanisms and their subsequent efficiency 
(Bekki 2009; McGarthy et al. 2008; Wetzel et al. 2013; 
Bekki 2014; Gen, Pop & Bahcall 2014; Bahe & McGarthy 

2015) . 

The uncertainty concerning environmental processing 
motivates a study in which we simulate in detail the funda¬ 
mental tidal and hydrodynamical mechanisms acting on a 
dwarf satellite. Specifically, we consider the apparent inef¬ 
ficiency with which 10^ M© satellites are environmentally- 
quenched (Wheeler et al. 2014, and references therein). By 
comparing the observed quenched fraction of these systems 
(Geha et al. 2012) with their associated fraction residing 
as subhalos in a more massive host (as deduced with a 
mock catalog constructed from the Millenium II Simula¬ 
tion), Wheeler et al. find that up to only 30 percent are 
quenched by their environment (most commonly a 
M© host). Using simple analytical models, they hypothesize 
that the corresponding quenching mechanism must act only 
on long timescales (e.g. ^ 7.5 Gyr). 

This inefficiency contrasts with the distinct transi¬ 
tion from blue irregulars to red early-types acknowledged 
amongst the smallest of the (10®“^ M©) dwarf galaxies when 
crossing the virial radius of the Local Group’s largest galax¬ 
ies (McGonnachie 2012). While ram pressure interactions 
with the Galactic hot halo and cosmic reionisation have 
been demonstrated to reproduce their current paucity of 
detectable gas (Mayer et al. 2001), the tidal stirring mech¬ 
anism (the repeated tidal shock/heating of a satellites stel¬ 
lar disc) reveals the efficiency with which groups environ¬ 
ments can transform very faint dwarf galaxies (Mayer et al. 
2001; Kazantzidis et al. 2011) and even Galactic-mass sys¬ 
tems (Bekki & Gouch 2011; Villalobos et al. 2012). Group 


mechanisms also appear to be responsible for the dramatic 
number growth, since at least z = 0.5, of the quenched 
and early-type systems known as Lenticular (SO) galaxies 
(Dressier et al. 1997; Willman et al. 2009; Just et al. 2010). 

In this analysis, we adopt methods demonstrated in our 
previous studies (Bekki & Gouch 2011; Yozin & Bekki 2014; 
Bekki 2014) and the results of previous cosmological studies 
that predict the infall epochs and group orbits for this stel¬ 
lar mass regime (e.g. De Lucia et al. 2012; Villalobos et al. 
2012), and devise an idealised model of this 10^ M© satel¬ 
lite and a 10^^ M© host system to characterise its long-term 
evolution. 

Other studies have considered the tidal trans¬ 
formation of dwarves within groups/clusters 
(Smith, Davies V Nelson 2010; Kazantzidis et al. 2011; 
Villalobos et al. 2012; Bialas et al. 2015), and the star 
formation histories of Galactic-mass satellites in groups 
(Feldmann, Garollo V Mayer 2011). Besides our novel 
choice of satellite and group masses, we distinguish 
this work with i) an estimation of the SF history and 
H I-deficiency using our chemodynamical model of the 
interstellar medium (ISM); ii) an interpretation of our 
results in terms of the stochasticity/evolutionary scatter in¬ 
troduced by infall inclinations/satellite interactions/group 
mass assembly, and hi) complementary simulations of 
hydrodynamical (ram pressure) interactions between the 
group medium and satellite. 

Section 2 describes the group and satellite models, while 
in Section 3, we illustrate the results of a suite of the 
tidal/hydrodynamical simulations. A summary and inter¬ 
pretation of these results in the context of previous obser¬ 
vational/numerical studies is given in Section 4. 


2 METHOD 

We model satellite galaxy infall into a 10^^ M© halo to assess 
the efficiency/mechanisms by which it is transformed and 
quenched. This involves placing an idealised dwarf galaxy 
(pristine bulgeless disc) at or beyond the virial radius (rvir) 
of a static group potential, with a trajectory prescribed from 
previous studies on the most cosmologically-common infall 
orbit. This evolution is assessed from a lookback time of 
8 to 10 Gyr, which represents the median epoch at which 
10® M© satellites resident in group-mass hosts (10^^ M©) 
at z = 0 were first accreted from isolation (McGee et al. 
2009; De Lucia et al. 2012). A schematic of this scenario is 
illustrated in Fig. 1. 

This method permits us to disentangle the primary 
mechanisms presently believed to govern dwarf satellite evo¬ 
lution, namely group tides, harassment, ram pressure and 
secular evolution. Grucially, we utilise satellite models with 
a 10^“^-fold improvement on baryonic resolution with re¬ 
spect to the cosmological-scale simulations whose results 
motivated this study. However, our method has disadvan¬ 
tages insofar as the group model neither follows the mass 
growth (by factor 2) expected from hierarchical growth over 
this timescale (Fakhouri, Ma V Boylan-Kolchin 2010), nor 
allows for the potential merging and/or resonant stripping 
of jointly accreted satellites (Gnedin 2003; d’Onghia et al. 
2009). We must further assume the instantaneous stripping 
of the satellites’ extended gas reservoir upon first crossing 
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Figure 2. Luminosity function for group satellites in the 
Tidal/ 2 ; = 1/Interactions simulations, shown in terms of 
binned luminosity (grey bars) and line-of-best-fit (red dashed 
line), compared with the equivalent data from SDSS DR4 
(Yang, Mo & van den Bosch 2008, ; Y08) for groups in the mass 
range 12.9<logio[M]<13.2, where open and closed circles denote 
satellite and central galaxies respectively 


Figure 1. Schematic of group models at 2 ; = 1 and 2 ; = 0 (left 
and right panels respectively). Black lines represent the satellite 
orbit in each case, with black circles denoting initial locations; 
red circles respresent the initial locations of other group satellites 
(size linearly correlated with satellite mass). In each case, the 
group halo density is conveyed in a logarithmic scale with the 
grey shaded region (darker=higher density). 


Table 1. Summary of group model and orbital parameters 


Parameter 

Value at ^ = 1 

Value at 2 ; = 

Virial mass (M©) 

10^3 

10^3 

Virial radius (kpc) 

329 

556 

Cnfw 

4.87 

9.74 

Apocentre (rvir) 

1.2 

1.0 

Pericentre (rvir) 

0.24 

0.2 

Pericentre Velocity (kms“^) 

650-700 

550-600 


the group ivir, such that this halo can be excluded from our 
modelling. This choice is motivated by the explicit hydrody¬ 
namic modelling of this process in galaxy-scale simulations 
(Kawata & Mulchaey 2008) and other studies which suggest 
the relatively dffuse intragroup medium (IGM) can even re¬ 
move a satellites’ cold gas disc (Bahe & McCarthy 2015). 

In the following sections, the group/galaxy models and 
simulations are described in further detail; we discuss our 
modelling choices in Section 4.2. 


2.1 Group Model 

Our group model is built on the assumption of a virialized 
group environment in which most substructure is accreted 
into a spherically symmetric background halo. For simplic¬ 
ity, both the dark halo and IGM are modelled with a NFW 
density profile (Navarro, Frenk & White 1996), which sta¬ 
tistical analyses of clusters suggest provides an acceptable 
mass distribution (van der Marel et al. 2000): 

“ r/rg,oup(l + r/rg.oup)2 ’ 

where r is the group-centric radius and rgroup = rvir/cNFW 
where cnfw is the concentration factor. 


To overcome the idealistic nature of our static group 
model, we adopt two representations of a group at distinct 
epochs of 2 ; = 1 and z = 0, whose combined results en¬ 
compass the characteristic evolution of a dwarf satellite ac¬ 
creted since z = 1. We fix the group mass to 10^^ M© at 
both epochs, and extract CNFw/rvir from the similar models 
of Villalobos et al. (2012); the adopted properties are sum¬ 
marised in Table 1. 

For the subset of simulations in which we explicitly add 
other group galaxies, we adopt the method introduced in 
Bekki V Gouch (2011) in which they are modelled as col¬ 
lisionless point masses, orbiting freely within the group, 
with an interaction softening length fixed to half each satel¬ 
lites’ rvir. The total mass of galaxies Mgr is established ac¬ 
cording to the total luminosity Lgr, with the relation of 
Marinoni V Hudson (2002): 


Ma 


= 350 


Ma^ 


5 X lO^^M© 


( 2 ) 


Monte-Garlo sampling is adopted to fit satellites, limited 
to a luminosity range 0.01 to 2.5L* and with an as¬ 
sumed mass-to-light ratio of 40, to the Schechter function 
with a faint-end slope -1.07 consistent with groups as es¬ 
tablished by Yang, Mo V van den Bosch (2008). Thus, our 
10^^ Md group hosts 52 satellites, with an initial spatial 
distribution prescribed with a NFW prohle and a con¬ 
centration of 3.0 to match the galaxy distributions re¬ 
vealed in K-band studies (Lin, Mohr V Stanford 2004). Fig. 
2 demonstrates good agreement between the correspond¬ 
ing linear ht to our LF with that obtained for groups in 
the mass range 12.9<logio[M]<13.2 from the SDSS DR4 
(Yang, Mo V van den Bosch 2008). The corresponding ve¬ 
locity dispersion of these satellites, a proxy for the inhuence 
yielded by satellite-satellite tidal encounters, is consistent 
with observations at between 125 and 150 kms“^. 


2.2 Orbits 

We derive the orbits from those exhibited by a ~ 10® satellite 
in the pure tidal/cohisionless simulations of Villalobos et al. 
(2012), who incorporated a live group halo to model its ca¬ 
pacity to retard the satellite via dynamical friction. A sat- 
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Time since infall (Gyr) 



Galaxy-centric radius (kpc) 


Figure 4. Rotation curve for the disc model (red solid line), as 
inferred from a/ GM (< r)/r, with separate curve for the stellar 
(red dashed line) component. The model is compared to recent 
3rd epoch HST data for the analogous LMC (black error bars; 
van der Marel & Kallivayalil 2013). 


Figure 3. A comparison of the orbits for a 10® Mq satellite in a 
10i3 group (red solid line), where the grey dashed line represents 
the host virial radius. Left panels show the simulated orbit for 
this scenario from Villalobos et al. (2012) in which the satellite 
is accreted with the cosmologically most common trajectory and 
is subject to dynamical friction by a live host halo; right panels 
shows the equivalent orbit (i.e. with similar mean pericentre radii 
and orbital period) used in this work. Top and bottom panels 
show the z = 1 and z = 0 scenarios respectively. 

isfactory live model of this halo requires significant reso¬ 
lution to avoid discreteness effects (a resource we choose 
instead to dedicate to the satellite model); given that an¬ 
alytical methods are not sufficiently accurate (Jiang et al. 
2008) and that this drag is nominally weak for dwarf satel¬ 
lites (Smith, Davies & Nelson 2010), we exclude this friction 
in our model. 

Fig. 3 conveys how we estabilish from Villalobos et al. 
an elliptical orbit at epochs z = 1 and z = 0 from the mean 
pericentres and apocentres across an 8-10 Gyr interval since 
first crossing the group rvir- In each case, the apocentre-to- 
pericentre ratio is consistent with the median among large- 
scale simulations (i.e. 4-5 Ghigna et al. 1998). Qualitatively, 
it can be seen that each orbit differs significantly in its peri¬ 
centre radii and orbital period; given that each is smaller 
in the z = 1 case, we expect its satellite evolution to be 
most dramatic. In each case, we deem the orbit sufficiently 
far from the group centre not to warrant the explicit mod¬ 
elling of a massive central galaxy often found lying at the X- 
ray temperature peak among groups (Zabludoff Mulchaey 
1998), and responsible for potentially efficient tidal shocks 
(Mayer et al. 2001). 

2.3 Disc Model 

Table 2 summarises the initial conditions of our satellite, a 
gas-rich bulgeless disc, which is comprised of a dark halo 
and stellar/gaseous disks. All galaxy models are comprised 
of 1.5x10® particles, divided between halo, stellar and gas 
components with 10®, 3x 10® and 2 xlO® particles respec¬ 
tively. 

For the stellar disc mass Md = 10^ M©, we find 
a halo mass Mh — lOOMd from observed intermediate-^ 
dwarf galaxies and abundance matching within large-scale 


Table 2. Summary of model parameters 


Parameter Value 


DM particle Number 1x10® 

DM-to-stellar mass ratio (M^/M©) 100 

NFW concentration (cnfw) 7-12 

Virial radius (rvir/kpc) 120 


Stellar particle Number 3x10® 

Stellar Mass (Md/M©) 10® 

Stellar scalelength (rd/kpc) 1.1 

Stellar scaleheight (^d/kpc) 0.22 

Toomre parameter (Q) 2.0 

Stellar metallicity ([Fe/H]) -0.9 

Metallicity gradient (d[Fe/H]/dr/dex/kpc) -0.04 

SNe expansion timescale (tsN/yr) 10®-® 


Gas particle Number 2x10® 

Gas mass fraction (fgas = Mg/Md) 1.0 

Gas scalelength (rg/rd) 2.6 

Gas scaleheight (zg/zd) 2.0 

SF Threshold Density {pth/cin~^) 10® 


simulations (Miller Bregman 2013; Munshi et al. 2013; 
Behroozi, Wechsler & Gonroy 2013), but note that the lat¬ 
ter may produce too much substructure, in which case our 
Mh may be a slight overestimate. We assume a top hat col¬ 
lapse model in which the mean halo density enclosed within 
rvir reflects the critical density at formation, and therefore 
use a scaling relation of the form rvir cx: M^® normalised by 
properties of the Galaxy (Klypin, Zhao Somerville 2002), 
from which we establish a satellite rvir of 120 kpc. We utilise 
a NFW density profile for the halo; the single distinction 
between our satellite at z = 1 and z = 0 lies in its initial 
concentration cnfw, which we constrain from mass-redshift- 
dependent relations of Munoz-Guartas et al. (2011), and es¬ 
tablish an evolution of cnfw since z = 1 from 7 to 12. 

An exponential disc morphology is adopted for the 
satellite, (truncated at 5 rd) with a density dehned for a 
cylindrical radius r and height z: 

/5(r,«) (xexp(-;^)sech^(-^), (3) 

rd 0.2rd 

where rd is a stellar scalelength derived from size-mass scal¬ 
ing relations (Ichikawa, Kajisawa <V Akhlaghi 2012) which 
exhibit a power law (index a = 0.15). This differs 
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from simple self-similar models of galaxy formation (e.g. 
Mo, Mao & White 1998) which do not account for mass- 
dependent processes such as outflow. Disc kinematics are 
constrained on the basis that the typical late-type disc 
with this Md will not host a secular bar instability (i.e. 
Melvin et ah 2014). Initial velocity dispersions, assigned ac¬ 
cording to epicyclic theory (with the vertical dispersion at 
a given radius is half that of the associated radial veloc¬ 
ity dispersion), are thus chosen with a Toomre’s param¬ 
eter of Q = 2 to inhibit expedient axisymmetric collapse 
(Combes & Sanders 1981). 

We enforce a gas-to-stellar mass ratio of unity from 
semi-empirical scaling laws for the gas mass (as a function 
of z) based on indirect measures of H I and H 2 content 
(e.g. the SF rate), which compares well with the high ob¬ 
served gas fractions of dwarves (Popping, Behroozi & Pee 
2014, and references therein). We assume this mass is also 
distributed with the aforementioned exponential profile, al¬ 
beit with a gas scalelength rg = 2.6rd as established from 
the normalised galaxy sample of Kravtsov (2013), and an 
assumed scaleheight Zg = 2zd. 

Fig. 4 illustrates how the initial rotation curve 
of this model fits acceptably to that of the LMC 
(van der Marel & Kallivayalil 2013), an exemplary galaxy in 
this stellar mass regime. 

2.4 Numerical modelling 

We utilise an original parallelised chemodynamical code 
GRAPE-SPH (GRAvity PipE-Smoothed particle hydrody¬ 
namics; Bekki 2009). Gravitational dynamics are computed 
with a softening length for halo and stellar particles (? 7 h= 2.1 
kpc and ? 7 s =200 pc respectively) which is set by the mean 
particle separation at the halo and stellar half-mass radii 
respectively. The ISM is modelled with smoothed particle 
hydrodynamics (SPH), with a minimum smoothing length 
of ~30 pc, and can cool radiatively (see Bekki 2014, for more 
details). This depends on an initial galaxy-wide mean metal- 
licity of [Ee/i7] = —0.9 prescribed from mass-metallicity re¬ 
lations (Kirby et al. 2013), and a fixed in-plane radial metal- 
licity gradient of 0.04 dex kpc“^. 

Star formation (SE) occurs via the conversion of gas 
particles to new stellar particles; this process is limited to 
collapsing Jeans-unstable regions with volumetric density 
exceeding a threshold pth, and relies on a probability pa¬ 
rameter, PsF, introduced in Bekki (2013): 

PsF = 1 - (4) 

where At is the timestep width for a given gas particle. If 
PsF exceeds a randomly generated number R with 0 ^ R ^ 
1, the conversion proceeds, based on the assumption that 
the timescale for star formation is shorter for higher density 
gas particles. This condition softens the numerical limitation 
that the giant molecular clouds which host SE lie below the 
resolution of gas particles. 

Stellar feedback is implemented with a sub-grid Super¬ 
nova (SN) model (see Bekki 2014, for more details) in which 
each subgrid SN type II expels a canonical 10^^ ergs, with 90 
percent delegated to thermal UV emission which is conferred 
upon neighbouring gas particles over a prescribed adiabatic 
expansion timescale tsn- The rate at which SNe occur is gov¬ 
erned by an initial mass function (IME), for which we adopt 



logpGas/(M0Pc')] 


Figure 5. Global (averaged) star formation law (H I-hH2 sur¬ 
face mass density vs. SFR surface density) for simulation Isolated 
(red circles), sampled at 0.25 Gyr intervals), compared with the 
Schmidt-Kennicutt law (power law index 1.4; black dashed line) 
for nearby normals and starburst spirals (Kennicutt 1998), and 
7 nearby spirals at 750 pc resolution (grey contours) reproduced 
from Bigiel et al. (2008). 

the standard Salpeter power law (a = —2.35) normalised to 
within the range 0.1 M© and 100 M©. 

This ISM model is simple in terms of its IME and the 
lack of an explicit modelling of H 2 -formation/outflow. How¬ 
ever, in Eig. 5 we illustrate how our choice of free parameters 
Pth and TsN (1 cm“^ and 10® yr respectively) can sufficiently 
replicate the observed SE law of dwarf galaxies Bigiel et al. 
(2008), who adopt far UV and 24muin emission to trace 
the inefficient (relative to the Kennicutt relation) SE in 7 
nearby spirals. We also ameliorate the concern that weak 
feedback schemes permit the existence of unphysical dark 
matter cusps, with the finding of a flat density profile within 
1-2 percent of the satellite rvu (e.g. Brook et al. 2012). 

Eor the subset of simulations in which we exclusively 
model the interactions between the satellite ISM and group 
IGM (i.e. no group potential), we adopt a method intro¬ 
duced in Bekki (2014) in which the self-gravitating galaxy 
model is placed within a cubic lattice of SPH particles (side 
length 25rd) that represents the IGM local to the satellite. 
The galaxy model is fixed, such that variations in the relative 
velocity and density of its surrounding IGM are applied by 
modifying at each time step the properties of the IGM SPH 
particles, according to the satellite’s orbit. Therefore, each 
IGM particle is assigned a density of ficMp('r), where r is the 
group-centric radius, p is the same NEW profile adopted to 
describe the group halo potential, and ficM is the mass quo¬ 
tient of the IGM with respect to the group halo. A fixed gas 
temperature of Ticm = 5.6 x 10® K is obtained for a 10^® 
M© group from Matsumoto V Tsuru (2000). 

2.5 Simulations 

This study analyses the results of 37 self-consistent simula¬ 
tions, whose details are summarised in Table 3. As a baseline 
case, we analyse the tidal influence of the group spherical po¬ 
tential z — 1 and z = 0 (labelled Tidal), with respect to 
the pure secular evolution of an identical galaxy (Isolated). 

Simulations labelled Ram refer to those in which a 
self-gravitating galaxy model interacts only hydrodynam- 
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Table 3. Summary of simulations 


Simulation 

Isolated 

Tidal/z = 1 

Tidal/^ = 0 
Ram/ 2 ; = 1 

Ram/z = 1/Overdense 
Tidal/ 2 ; = 1/Orientations 

Tidal/ 2 ; = 1/Harassment 


Description 

Isolated model (not orbiting 
within group potential) 
Satellite under tidal influence 
from the group potential only; 
orbit and group properties at 
z = 1 

As above, but with orbit and 
group properties at 2 ; = 0 
Satellite under Ram Pressure 
only; properties at 2 ; = 1 
As above, but with IGM den¬ 
sity increased by factor 10 
Tidal/ 2 ; = 1, but with 16 vari¬ 
ations on the initial orientation 
w.r.t orbital plane 
TidaR = 1, but with other 
group satellites included, 16 
variations on initial position 


ically with the IGM. We conduct two simulations to ad¬ 
dress the variations in the local density experienced by in¬ 
falling satellites. The first, Ram/z = 1, adopts the z = 1 
group properties and a ficM of 0.1 consistent with the mean 
group ICM mass inferred from X-ray observations (Lin et al. 
2012). In this case, our adopted ICM density profile com¬ 
pares sufficiently well with those which fit the distribution 
of galaxies, as inferred for example from X-ray temperature- 
satellite velocity dispersion relations (Zabludoff & Mulchaey 
1998). A second simulation, Ram/z = 1/Overdensity, 
adopts ficM = 1 to replicate regions of live hydrodynamic 
models of the cluster medium local to an infalling satel¬ 
lite (e.g. Tonnesen & Bryan 2008), which exhibit a mag¬ 
nitude enhancement in p with respect to the spherically- 
averaged p at that radius, and are the preferential regions in 
which quenched satellites can be located (Bahe & McCarthy 
2015). 

Simulation Tidal/z = 1/Harassment denotes the com¬ 
bined results from 16 simulations in which we model the 
harassment of our infalling disc satellite galaxy by other 
group satellites (modelled as point massses), in addition to 
the underlying group tidal field. To appraise the diversity 
and stochastic nature of harassment, each simulation differs 
in the prescribed initial location of the infalling satellite. 
In each case, this location is randomly selected from the 
locus defined by the group rvir- Although simplistic in its 
approach, the mean and la variation in these simulations’ 
evolution obtained from this Monte-Carlo method provides 
a useful comparison with our baseline case in determining 
the systemic effect of harassment. 

Simulation Tidal/z = 1/Orientations denotes the com¬ 
bined results from 16 simulations in which the initial incli¬ 
nation of the satellite with respect to its orbital plane is 
assigned randomly. Again, these results are compared with 
our baseline cases in which an intermediate infall inclina¬ 
tion is adopted, comprising a -45 degree angle between the 
angular momentum vector of the disc and the group z-axis, 
a -30 degree azimuthal angle between group x-axis and the 
projection of the disc angular momentum vector onto the 
group x-y plane (corresponding to a prograde orbit). 


2.6 Diagnostics 

Simulation data is stored and reduced at 70 Myr increments. 
The following diagnostics are utilised to compare the evo¬ 
lution of each simulation in terms of transformation and 
quenching/gas deficiency. 


2.6.1 Morphology 

The stellar/DM/gas mass dynamically associated with the 
satellite is established with an iterative procedure, wherein 
gravitationally-bound material discerned in one step is used 
to build a more refined spherical mass distribution of the 
galaxy in the following step, until the centre of bound mass 
is converged upon. 

Stellar magnitudes are approximated in this 
work with synthesis models for the evolution of stel¬ 
lar mass to light ratios (Bruzual & Chariot 2003; 
Portinari, Sommer-Larsen & Tantalo 2004). Assuming 
infall at a lookback time of 8-10 Gyr, an initial stellar age 
distribution is established from an estimated SF history 
up to 2 ; = 2.5, calculated using extrapolated SFR-stellar 
relations from Whitaker et al. (2012). We approximate the 
H I -mass as a factor 1.38 ^ of each gas particle mass. The 
attenuation by dust is approximated with the gas column 
density and an assumed gas-to-dust-extinction ratio of 
5x10^^ cm“^Av, consistent with bright dwarf/LMC-type 
galaxies (Dobashi et al. 2008). 

The surface brightness profile p(r) is composed from the 
corresponding attenuated luminosity, from which a bulge-to- 
disc mass ratio (B/T) is estimated with the excess light lying 
within 2 effective radii with respect to a fitted exponential 
disc profile, and normalised with the integrated flux of the 
disc. The fitted disc profile is estimated with an algorithm 
that progressively smoothes the radial light profile until the 
underlying disc (which we find to exist in all cases except 
where the disc is completely destroyed) can be distinguished 
from the bulge and/or (positive or negative) truncation com¬ 
ponents of p{r) by lying coincident with a global minimum 
of dV/dr^ 

A B/T^O.5 is often associated with an early-type mor¬ 
phology (Kormendy & Kennicutt Jr. 2004; Willman et al. 
2009), although recent studies have highlighted an inherent 
Haw in adopting simple fits to /i(r) insofar as it incorpo¬ 
rates features such as bars/ovals/lenses in the excess light. 
In the case of early-type systems such as dSO/SOs in which 
these features are hypothesized to be critical to their for¬ 
mation, this can lead to a significant overestimate of B/T 
(Laurikainen et al. 2009). 

For this reason, we also compare our satellite’s stellar 
kinematics with those of observed E and/or SOs, using a 
specific angular momentum Ar computed in the manner of 
Emsellem et al. (2007), who adopt Ar to distinguish galaxies 
as fast/slow rotators: 

A„ - . (5) 

where F^, Ri, Vi and ai are the total flux, radius, mean 
stellar velocity and velocity dispersion respectively, at the 
zth radial bin up to the effective radius (re) 

A quantitative demarcation between fast- and slow- 
rotators can be asserted by the relation A = 0.31 y/e found 
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Figure 6. Snapshots from simulations Tidal/^ = 1, Tidal/ 2 ; = 0, Ram/ 2 ; = 1 and Ram/ 2 ; = 1/Overdense. The top row illustrates the 
B-band surface brightness of Tidal/ 2 ; = 1, viewed edge-on, to convey the thickening of the disc (with time since first infall labelled above, 
and the axial ratio c/a also provided). For all rows below, the underlying colour scale corresponds to the surface gas mass density, overlaid 
with an ellipse (white dotted line) fit to the 28 mag arcsec“^ isophote of the B-band surface brightness. Tidal simulations snapshots are 
provided with the instantaneous value of r^ar? An and defni; the bottom two rows show Ram simulations with the instantaneous defni- 


among 260 early-type galaxies in the ATLASproject 
(Emsellem et al. 2011), where the the apparent stellar ellip- 
ticity £ = 1 —0.5c/a, and the axial ratio c/a is computed here 
from the eigendecomposition of the angular inertia tensor for 
all bound stars with r < 2re. Approximately, Ar ^ 0.2 cor¬ 
responds to a bright early-type system (E), while SOs and 
fainter ellipticals can exhibit significant rotational support; 
this distinction can provide insight as to their formation, 
featuring less overlap between types than other metrics e.g. 
V/cr. 

The formation/survival timescales of bisymmetric 
structures that are associated with transformation mecha¬ 
nisms are conveyed with two metrics. The bar radius, nor¬ 
malised by the stellar scalelength (rbar/ra), is estimated via 
an empirical method (see Yozin & Bekki 2014, and refer¬ 
ences therein); briefly, we detect a bar if a succession of 
ellipses fit to isophotes of the deprojected, face-on surface 
brightness distribution (incrementally, from 21 to 28 mag 
arcsec“^) deviate in phase by less than 10 degrees, where 
the outer most ellipse has an axis ratio no larger than 0.4. 

The disc-wide presence of spiral/tidal arms is quantified 
with An, the mean amplitude of the second Eourier mode 
in concentric shells spanning radii from the bar radius up to 
2 effective radii; where spiral arms typically register in the 
range 0.15<Aii<0.6 (Rix & Zaritsky 1995). 


2.6.2 Star formation 

We express the logarithmic SE rate (SER) in terms of an 
offset (ASER) from that expected for a galaxy of its in¬ 
stantaneous stellar mass, for which we adopt the linear rela¬ 
tion devised from star-forming galaxies in the SDSS catalog 
by Kauffmann et al. (2003b). By considering the colour bi¬ 
modality expressed by the general galaxy population, a pas¬ 
sive classification for a galaxy can be determined if ASER<- 
1 (i.e. SER reduced by a factor 10), with the red sequence 
lying at approximately -2 (see also Bluck et al. 2014). A sim¬ 
ilar result is obtained if we assume that a galaxy can be 
classified as quenched if its specific SER (sSER) ^ 10“^^ 
yr“^, as implied from a bimodality in the sSERs of central 
galaxies in the SDSS DR 7 (Wetzel et al. 2013). 

Eurthermore, the visual signature of SE can manifest in 
strong Ha emission, a tracer of massive star formation which 
we can approximate with the relation (Kennicutt 1998): 

L{Ha) = SFR{MQyr~^) x 1.26 x 10^\rgs~^ (6) 


where the SER is estimated from new stellar particles with 
age less than 3x10^ yr (the approximate age of an ~8 Ma 
star). 
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Figure 7. Diagnostic of the evolution of simulations Tidal/ 2 ; = 1 (solid red line), Tidal/ 2 ; = 0 (dotted red line) and Isolated (dashed 
black line). (Clockwise from top left panel) Star formation rate (logarithmic) expressed as an offset from the expected SF for its stellar 
mass (from Kauffmann et al. 2003b), where a passive or quenched classification is taken to lie at an ASFR of -1 or -2 respectively; 
Bulge-to-total luminosity ratio (B/T); stellar bar radius r^ar? normalised by the stellar scalelength r^; amplitude of the second Fourier 
mode An averaged in radial bins up to the tidal radius, where the signature of spirals lie in the range 0.15 (weak arms) to 0.6 (strong); 
specific angular momentum of all stars within the effective disc radius Ar; H I-deficiency of the galaxy, describing the H I-mass within 
its optical diameter with respect to that expected for its morphological type as established from B/T. 


2.6.3 Gas mass associated with the satellite 

We quantify the simulated loss of gas with a defei parameter 
(Solanes, Giovanelli & Haynes 1996): 

defHi = a + 61ogio(/iDopt) - logio(Mjf (7) 

where h is the Hubble parameter (~0.7) and Dopt is the 
linear optical diameter, which we assume equivalent to the 
galacto-centric radii at which the deprojected azimuthally- 
averaged surface brightness falls under 25 mag arcsec”^. 
Mhi is the total mass of gravitationally bound gas in our 
simulated galaxy, and a and b are statistically-inferred co¬ 
efficients which vary according to morphological type. At 
each time step of our simulations, we classify the type with 
B/T and linearly interpolate between the appropriate coef¬ 
ficients; this, combined with a D 25 that evolves according to 
our adopted stellar synthesis models, yields a defni param¬ 
eter comparable with observations. 


3 RESULTS 

3.1 Dependence of the satellite evolution on 
group tides since infall at 2 ; = 1 

Fig. 6 shows the deprojected face/edge-on distribution of 
stars (in /ie) and gas from our baseline Tidal and Isolated 
simulations, in a series of snapshots spanning 10 Gyr since 
their first crossing of the group virial radius. It illustrates 
how their first pericentre following infall excites the forma¬ 
tion of tidal arms and a bar; the associated tidal heating 
thickens the disc (resulting in boxy isophotes), while the in¬ 
corporation of an ISM, funnelled to the centre by the bar 
and tidal torquing, contributes to bulge growth and bar 


dissipation. Neither tidal stripping nor secular/induced SF 
can remove the satellites’ gas completely within this 10 Gyr 
timescale. 


3.1.1 Star formation history 

In the absence of gas accretion, the SF rate among our base¬ 
line simulations is shown in Fig. 7 to follow an exponential 
decline in SFR. This is interrupted only by tidally-triggered 
SF at each pericentre passage, with an enhancement fac¬ 
tor of up to 5 consistent with previous numerical studies 
(Di Matteo et al. 2008), and attributed to gas infall from 
the H I-rich disc. 

The ASFR metric implies the mean timescale of 8 Gyr 
since the accretion of a luminous dwarf galaxy is sufficient 
for its classification as passive, with a stellar mass growth in 
the intervening time of only 20 percent. From SFR-mass re¬ 
lations (Whitaker et al. 2012), a characteristic SFR of order 
1 M 0 yr“^ betwen 2 ; = 0 — 1 for a 10^ M© galaxy would in¬ 
dicate an analogous isolated galaxy could grow by a factor 
of 1-5, the majority of which fuelled by a cold gas accre¬ 
tion mode which dominates in this stellar mass regime for 
all 2 ; (Dekel et al. 2009). In conservative analytical models 
approximating the outflow strength and recycled fraction of 
gas, the SFR due to fresh accreted gas can be as much as 
50 to 90 percent of the total SFR within a gas consumption 
timescale, for those galaxies hosted by low mass halos in 
which outflow is relatively efficient (Sanchez Almeida et al. 
2015). The deficit between the initial SFR of our star form¬ 
ing simulations with respect to ASFR=0 can therefore be 
explained by the exclusion of accretion and a gaseous halo 
from our code, but the subsequent evolution will depend 
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Figure 8. (Left panels) Surface mass density of new stars (top) 
and Toomre stability parameter for the gas disc (Q) as a function 
of disc cylindrical radius, at 8 Gyr after infall, for Tidal^ = 1 (red 
line/empty circles) and Isolated (black line/filled circle) simula¬ 
tions; (Right panels) Ha emission density, in logarithmic scale, 
for the same simulations/time. 

sensitively on the timescale in which this accretion is shut 
off, which we here assume as instantaneous. 

An observable implication of this unreplenished gas re¬ 
serve, combined with the influence of group tides, is the 
outside-in truncation of SF (Fig. 8 ). Previously reproduced 
by Bekki & Couch (2011) in simulations of massive group 
spirals and attributed to the heating of the outer disc/gas 
infall, we extend this result here to dwarf satellites, and 
indeed find the radial extent of SF typically lies where 
the stellar/gas Q parameters (= cr/^/3.36GE, where a, k, 
G and S are the radial velocity disperion, epicyclic fre¬ 
quency, gravitational constant and surface mass density re¬ 
spectively) exceed the canonical range of 1.5-2 required for 
axisymmetric stability. A similar phenomenon was obtained 
by Bekki, Couch & Shioya ( 2002 ), who also attributed the 
growth in Q and concurrent spiral arm fading with ceased 
halo gas accretion. 

The ongoing presence of Ha emission beyond 8 Gyr ago 
infall means these satellites (under exclusively tidal influ¬ 
ence) would not meet the definition of a quenched system 
(e.g Geha et al. 2012 ) nor that of a classical SO; their blue- 
cores are however comparable with a sample of dEs/dSOs 
observed in recently accreted subgroups of the Ursa Major 
cluster (Pak et al. 2014) for which no evidence of significant 
hydrodynamical influence is found. 

3.1.2 Gas and halo loss 

An initial defni of zero is consistent with our satellite com¬ 
mencing with the mean structural properties/fgas of its 
mass/type, with a subsequent increase with time corre¬ 
sponding to tidal stripping, outflow and consumption of 
gas via SF. Uncertainties in defni among observed galax¬ 
ies can be large (~0.25 dex; Pappalardo et al. 2012), mo¬ 
tivating Cortese, Catinella & ad S Heinis ( 2012 ) to adopt 
defni>0.5 to denote a deficient galaxy among their sample 


Figure 9. (Main panel) For simulations Tidal/ 2 ; = 1, Tidal/ 2 ; = 0 
and Isolated, each band conveys the gas mass bound to the galaxy 
(lower limit) and the gas mass, plus stellar mass formed from 
gas, bound to the galaxy (upper limit), all normalised by the 
initial gas mass, as a function of simulation time; (Inset panel) 
the normalised DM mass bound in Tidal/ 2 : = 1 (black solid line), 
and the normalised total mass within the tidal radius (red dashed 
line). 

of 322 nearby galaxies. Adopting this threshold, we find a 
timescale of 8 Gyr within the group environment sufficient 
to render an originally gas-rich galaxy H I-deficient, albeit 
not removed of gas entirely. 

Fig. 9 shows how this reflects a loss of at least 30-60 
percent of its original cold disc gas; by comparison, the dark 
matter halo is truncated such that only 5 percent of its orig¬ 
inal mass remains bound to the galaxy (Fig. 9); In practice, 
the dynamical mass evolution within the tidal radius of the 
satellite (which is dominated by baryons) only amounts to a 
largely monotonic reduction in time by 50 percent. This re¬ 
sult contrasts sharply with those of Wetzel et al. (2013), who 
propose a substantial stellar mass growth following group 
infall as part of a two-stage quenching process wherein the 
initial stripping of the satellite’s halo gas is inefficient. 

The gas mass loss exhibited in Tidal simulations in¬ 
cludes a secular component which we assume equal to that 
demonstrated by Isolated; this amounts therefore to about 
20 percent loss over 10 Gyr, divided between consumption 
by SF and outflows. As a sanity check, the corresponding 
gas depletion timescale of Isolated (~50 Gyr) is compared 
with that expected for its stellar mass, where for the latter 
we assume a typical Hl-to-stellar mass ratio for a 10^ M© 
galaxy of 1 to 5 (e.g Haynes et al. 2011) and a median SFR 
of 0.1 to 0.2 M 0 /yr (Kauffmann et al. 2003b). This yields an 
acceptable discrepancy by a factor 2-3, broadly attributed to 
our simple outflow model and the artifical closed-box nature 
of Isolated. 

3.1.3 Transformation mechanisms and efficiency 

Fig. 10 illustrates the light profile of Tidal/z = 18 Gyr after 
initial infall. In contrast with the single exponential profiles 
obtained in the collisionless but otherwise similar simula- 


MNRAS 000 , 1-16 (2013) 













10 C. Yozin and K. Bekki 



Figure 10. (Main panel) Surface brightness profile, as a func¬ 
tion of disc cylindrical radius, at 8 Gyr after infall for simulation 
Tidal/^ = 1. The best-fit line to the disc component is shown 
(black dashed line); the B/T ratio is computed from the excess 
light with respect to this line, lying within the tidal radius (rtidal)5 
(Inset panel) a comparison of /^-profiles for Tidal/z = 1 at 4 and 
8 Gyr after infall. 


tions of Villalobos et ah (2012), our satellites’ profiles com¬ 
prise a bulge component strongly linked with tidally-induced 
SF (Christlein V Zabludofit 2004), and an anti-truncation 
(up-turn) in the outer disc whose origin is less clear. We 
find tentative evidence for the break strength T (ratio of 
inner and outer disc scalelengths) to increase as a function 
of time resident in the group potential (inset panel). This 
contrasts with an analysis of the V-band /i(r) profiles of 280 
classical SOs by Maltby et al. (2015), a morphological type 
also preferentially formed in groups (Willman et al. 2009) 
but for which no environmental-dependence of T is found. 
Instead, Maltby et al. compare their sample with equivalent 
spiral galaxies and propose this up-turn is a manifestation 
of extended bulge light that dominates a fading disc. 

The quantitative growth of the central bulges of 
Tidal/z = 1/Tidal/z = 0 and Isolated is shown in Fig. 
7. In the absence of merger activity, these are not clas¬ 
sical bulges, as further indicated with line-of-sight veloc¬ 
ities of the disc which reveal a cylindrical profile in the 
nucleus. The slow secular growth exhibited by Isolated is 
consistent with the assertion by Kormendy V Kennicutt Jr. 
(2004) that pseudo-bulges can comprise up to 10 percent of 
the total disc mass. By contrast, the bulge growth of the 
Tidal simulations, whose B/T is consistent with the corre¬ 
sponding range, for this stellar mass, of 0<B/T<0.4 collated 
by Bluck et al. (2014), emerges principally as a series of large 
enhancements coincident with pericentre passages. 

The underlying mechanism of this bulge growth is re¬ 
vealed by the formation of a stellar bar, as indicated by 
fbar/rd- In both Tidal simulations, the first pericentre pas¬ 
sage following infall is sufficient to tidally-trigger this struc¬ 
ture to its largest observed radius. A subsequent enhance¬ 
ment in the SF rate and increase in B/T is thus consistent 
with gas infall being faciliated by tidal torques/the stellar 
bar. However, this simple observable correlation between 
bar and bulge breaks down with time because the growth 
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Figure 11. Evolution of specific angular momentum (within 
re) for simulations Tidal/z=l and Tidal/z=0 (red plus and circle 
symbols respectively) as a function of stellar ellipticity. Points 
are plotted at 0.5 Gyr increments; as an idealised bulgeless disc, 
each simulation commences in the vicinity of (e = 0.8, Ar = 0.8). 
Black bars illustrate the Icr distribution in each dimension for 300 
GALIFA galaxies classified by morphological type (reproduced 
from Falcon-Barroso, Lyubenova & van de Ven 2014). The grey 
arrow conveys the general evolution of our Tidal simulations; the 
dotted line is an effective demarcation between fast/slow rotating 
galaxies defined by Ar = (Emsellem et al. 2011). 


- Simulation ^ 
time 

Sb- 


Sa- 




so 




(Stoe °+ 


+ 

+ + 


^ o o. 


Sc 


.0.31 Ve 


Sd 


Tidal/z=1 + 
Tidal/z=0 o 

J_I_I_I_I_ 


of the bulge increasingly dominates the shape of isophotes 
in the inner disc and thus weakens the visual signature of 
the bar, whilst the concentration of SF in the bulge means 
B/T can also grow due to the relative fading of the bar and 
disc. A decline in rbar/ra is also to be expected after several 
pericentre passages due to the buckling of the bar by sig¬ 
nificant gas infall, similarly reported in simulations of Virgo 
clusters satellites by Mastropietro et al. (2005), and the pro¬ 
gressive weakenening of the bar through angular momentum 
exchange (Bournard, Combes V Semelin 2005). 

An important result however is that the simulated B/T 
is at no point consistent with the range 0.5^B/T associated 
with an early-type morphology (Kormendy V Kennicutt Jr. 
2004; Willman et al. 2009). This limited morphological evo¬ 
lution, wherein late-type structure is fundamentally main¬ 
tained, is also reflected in Ar. In expanding upon Fig. 7, 
which highlights the sharp loss in angular momentum upon 
each pericentric passage of the Tidal simulations (with re¬ 
spect to the stability of Isolated), Fig. 11 illustrates the evo¬ 
lution of Ar as a function of the apparent stellar ellipticity 
£. Their evolution towards lower e and Ar is consistent with 
the thickening/heating of the inner disc by tidal shocks upon 
pericentre passages and the triggered bar. The final kine¬ 
matics thus ultimately resemble those observed for Sa and 
SO morphological types, lying in agreement with our prior 
determination of its type with B/T. Gas accretion from fil¬ 
aments/other galaxies could in principle raise Ar if corotat¬ 
ing with the stellar disc, although we have assumed no such 
process occurs within our group enivronment. 

A key feature of these Tidal simulations that precludes 
their association with the transitory SO classification is that 
disc spiral structure beyond the bar proves to be quite ro¬ 
bust, being triggered by tides upon all four pericentre pas¬ 
sages of the Tidal/z = 1 simulation. The peak amplitudes 
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of these arms are short-lived (or order a Gyr) and weaker 
on successive orbits because i) the stellar disc is increasingly 
concentrated, and ii) SF that would otherwise be triggered 
in these outer disc regions is suppressed by the prior in¬ 
fall/stripping of gas and the absence of external gas accre¬ 
tion (e.g. Bekki, Couch & Shioya 2002). Nonetheless, there 
is only a ~Gyr window (in Tidal/z = 0) in which our Tidal 
satellites do not bear the signature of an arm (i.e. where 
An < 0.15), strongly suggesting they cannot in general ex¬ 
hibit the axisymmetries inherent to E/SOs. 

3.2 Dependence of the satellite evolution on 
orbital inclination 

Fig. 12 compares the baseline simulation Tidal/z = 1 with 
the combined results of Tidal/z = 1/Orientations, wherein 
we find the former corresponds closely to the mean evo¬ 
lution of the latter in the various metrics that describe the 
satellite evolution. Variations, in particular those of the trig¬ 
gered starbursts upon first infall, primary stem from how the 
efficiency with which tidal shocks can trigger bar instabili¬ 
ties/gas infall depend on the disc/orbital plane alignment. 
On long timescales, however, this dependence becomes am¬ 
biguous (e.g. due to stochastic events/resonances of the disc; 
Kazantzidis et al. 2011). Moreover, the stripping efficiency 
of the extended gas disc and DM halo (the latter being pre¬ 
scribed with no initial net angular momentum) are (to first 
order) invariant with respect to infall inclination. 

The clear exception lies in the growth of B/T which, 
although in part a manifestation of tidal heating (due to 
tidal shocks which prove weakly dependent on inclination), 
reflects also a sensitivity to the efficiency of bulge growth 
by the starburst upon first infall (Ghristlein & Zabludoff 
2004). Acting upon a near pristine stellar/gas disc, the first 
group tide interactions that can induce a strong bar/gas 
infall in a prograde edge-on/inclined scenario (adopted in 
Tidal/z = 1) are conversely inefficient in a pure face-on 
and/or retrograde scenario. However, a weak enhancement 
in BjT in these latter cases corresponds to an inefficient 
manner in which the associated mechanisms described in 
Section 3.1.3 can destroy/obscure the bar. 

We therefore identify a complex relationship between 
infall inclination and structural properties of a star-forming 
group satellite, which manifests in a broad range of H/T 
(since a mean infall lookback time of 8 Gyr) consistent with 
observations for this mass range (Bluck et al. 2014). How¬ 
ever, the transition to a quenched status over this timescale, 
reflected in ASFR and defni, remains well described by the 
reference Tidal models of Section 3.1. 


3.3 Dependence of the satellite evolution on 
satellite-satellite interactions 

Fig. 12 compares the baseline simulation Tidal/z = 1 with 
the combined results of Tidal/z = 1/Harassment. The prin¬ 
ciple remark from this comparison lies in the broad scat¬ 
ter revealed by the latter sample; the width of this scat¬ 
ter depends on the interaction softening length, although 
our adopted constraint is compatible with the typical min¬ 
imum separation between encounters (~50 kpc) found by 
Villalobos, De Lucia V Murante (2014), who similarly es¬ 


tablish that tidal disruption of a satellite on an infalling or¬ 
bit is, to first order, governed by interactions with the group 
halo rather than harassment by other group members. 

Satellite encounters at large orbital radii trigger an el¬ 
evated SF rate upon first infall compared with Tidal/z = 1 
(and an associated enhancement in defni also contributed 
to by tidal stripping), although the SF history largely re¬ 
mains defined by the secular exhaustion of the original gas 
disc. As noted in Section 3.1.3, the earliest triggers of gas in- 
fall/starbursts largely define the bulge mass, as subsequent 
tidal interactions are unable to trigger as large/strong a bar. 
These interactions are not sufficient to disrupt the disc di¬ 
rectly; instead, three-body interactions or those of unequal 
mass ratios can either send the dwarf satellite to a higher 
energy orbit (increased pericentre, period) where the tidal 
influence of the group halo is weaker, or towards a efficiently 
disruptive path to the group centre. Two such examples are 
illustrated in Fig. 12 with a highly divergent evolution in 
their respective B jT. 


3.4 Dependence of the satellite evolution on 

interactions with the Intra-Group Medium 

Fig. 13 compares the evolution in ASFR and fgas as exhib¬ 
ited by Isolated and Ram simulations. A key observation 
here lies in the enhancement in SF wth respect to Isolated 
following first infall, a phenomenon similarly reported by 
Bekki (2014) who attributed this star burst to the compres¬ 
sion of the satellites’ ISM at the IGM-ISM interface (see also 
Kronberger et al. 2008). In Ram/z = 1, this enhancement 
can be long-lived because the stripping of the satellites’ cold 
disc gas is inefficient in this scenario, with a final fgas of 60 
percent. 

The evolution of simulation Ram/z = 1/Overdensity 
indicates that a localised pocket of overdense IGM is suffi¬ 
cient to remove ~50 percent of the satellite’s cold gas disc 
within only its first orbit of the group. If combined with 
the gas loss exhibited in Tidal/z = 1 (i.e. combining simul¬ 
taneous tidal and ram pressure mechanisms acting on the 
satellite since first infall), this satellite can lose all its gas 
and therefore be effectively quenched within 4-5 Gyr of first 
passing rvir- 

This result is consistent with an apparent peak in the 
frequency of satellites being quenched by ram pressure strip¬ 
ping at z = 1, where this frequency is not simply a mono¬ 
tonic function of time (Bahe V McGarthy 2015). However, 
our result is largely qualitative insofar as our idealisation 
of the IGM distribtion of a typical group oX z — \ carries 
some uncertainities. The strong dependence of the quench¬ 
ing timescale on the local IGM density of the satellite is 
therefore the main result of this section. 

In an attempt to verify the gas loss simulated numer¬ 
ically in Ram/z = 1, we draw upon the analytical formu¬ 
lation for the satellite stripping radius (rgtr, at which pV^ 
exceeds the orthogonal restoring force of the galaxy disc) of 
Gunn V Gott (1972): 

py2>S(r)^, (8) 

where S(r) is the gas surface mass density at galacto-centric 
radius r, and 6<j)l6z is the gradient of the potential field 
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Figure 12. Diagnostic of ASFR, defni, and B/T (as in Fig. 7), shown for simulation Tidal/ 2 ; = 1 (red solid line), and the mean/lcr 
variation (black dotted line/grey shaded region) in the respective evolution of two sets of simulations run in accord with Monte-Carlo 
sampling, where the top row shows the results for 16 variants of Tidal/ 2 : = 1 where the initial orientation of the satellite (with respect to 
its orbital plane) is assigned randomly (Orientations), and the bottom row shows the results for 16 variants of Tidal/ 2 ; = 1 where other 
group galaxies are incorporated in the model as point masses and the primary satellite is assigned a random initial location on the locus 
defined by the group halo rvir (Harassment). In the B/T panels, the evolution of selected individial simulations from each set are shown, 
including a satellite accreted face-on in inclination, and two simulations in which harassment strongly alters the primary satellite’s orbit 
and thus its transformation process (red dotted lines). 


perpendicular to the disc plane. While idealised for a spi¬ 
ral galaxy with co-axial orbital velocity and angular mo¬ 
mentum, Roediger & Hensler (2006) demonstrate that this 
method can be extended to an arbitrary orientation such as 
that adopted here. We compute rgtr per time-step using the 
(DM, stellar and gas particle) mass distribution and local 
piCM of our simulated galaxy, and adopt the maximal 54>/6z 
to account for the warp of the disc with respect to our plane 
of best fit and other difficulties arising from this analytical 
method as discussed in Roediger & Briiggen (2005). By ap¬ 
plying this rstr to the initial structural definition of the gas 
disc, we compute an analytical fgas from the accumulated 
gas mass lying beyond rgtr, which we assume instantaneously 
removed. 

Illustrated in Fig. 13, this analytical fgas slightly overes¬ 
timates the gas loss simulated in Ram/z = 1. This partly re¬ 
flects the capacity for stripped gas to not actually exceed its 
escape speed, but more pertinently stems from an assump¬ 
tion in the analytical method for a maximal ram pressure 
~ where V is the satellite velocity Vgat- In practice, 
the ICM envelope is not stationary, and the galaxy velocity 
and relative ICM-galaxy velocity can be strongly correlated. 
A similar result was reported by Tonnesen & Bryan (2008) 
who establish, for a satellite pericentre velocity Vgat of ~500 
kms“^ (comparable to Ram/z = 1) in a cluster-scale envi¬ 
ronment, the ICM ahead of a satellite’s trajectory can be 
travelling at a rate of as much as 0.5Vgat- 


4 DISCUSSION AND CONCLUSION 

In Section 2, we developed a representative model of a gas- 
rich late-type dwarf galaxy (stellar mass ~ 10^ M©), using 


scaling relations and the mean properties of isolated galaxies 
in this mass regime, to match their observed mass distribu¬ 
tions and SFR as a function of gas surface density. In Sec¬ 
tion 3, we establish the typical chemodynamical evolution of 
this satellite following first infall to a (10^^ M©) group envi¬ 
ronment in a suite of simulations. We address the idealistic 
nature of the group model (where significant mass assembly 
is not modelled) by characterising the satellites’ evolution 
as bounded by two group models/infall orbits at 2 ; = 1 and 
z = 0. 

We first establish, in a pure tidal model of a group (Sec¬ 
tion 3.1), that the mean time in which a low- 2 ; satellite of this 
mass has been resident with a group (~8 Gyr) is sufficient to 
suppress SF to a passive (and H I-poor) status, if assuming 
rapid/instantaneous removal of external gas sources (Fig. 
6). This galaxy is not quenched, insofar as SF is ongoing 
(but centrally concentrated); tidal heating plays an active 
role although the morphological transition to an early type 
is not complete within a 8-10 Gyr timescale. 

In sections 3.2-3, we take a Monte-Garlo approach to 
assert the sensitivity of this evolution to a variety of orbital 
inclinations and satellite harassment (Fig. 11), and confirm 
that the above result is indeed sufficiently characteristic of 
a satellites’ typical evolution. There exists however a broad 
scatter, in particular for those simulations in which other 
group members are included, with the implication that the 
characteristic quenching/transformation can be mitigated 
for a significant proportion of infalling satellites whose orbits 
acquire significant energy. 

In section 3.4, we simulate numerically the hydrody- 
namical ISM-IGM interaction (ram pressure stripping) with 
results that are consistent with a classic analytical model 
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Figure 13. (Top panel) ASFR for simulations Isolated, Ram/z = 
1 and Ram/^ = 1/Overdensity; (bottom) fgas (=Mhi/Mhi, initial) 
for simulations Isolated, Ram/z = 1, Ram/ 2 ; = 1/Overdensity 
and Tidal/^ = 1. Also shown is an analytical estimate of fgas 
(black dotted line), using the approximation of rstr for the z = 1 
orbit and assuming an exponential cold gas disc. 


(Fig. 13), albeit with deviations that can be explained in 
terms of a correlation of satellite and local ICM velocities. 
We find in general an enhanced starburst upon first infall, 
followed by a suppression whose rate depends on the local 
ICM density; if matching the density to that preferentially 
found around quenched satellites of previous numerical stud¬ 
ies, we find a significant gas loss imposed by ram pressure 
stripping, with a corresponding quenching of SF within 4-5 
Gyr of first infall at z = 1. 

To conclude, our simulations are consistent with ob¬ 
servational evidence for the inefficiency with which a 10^ 
M© satellite is quenched in a (10^^ M©) group environment 
(Wheeler et al. 2014), if most such satellites maintain their 
relatively slow orbital decay following satellite harassment, 
and are not resident within a significant local overdensity of 
the intra-group medium. Our simulations are however ide¬ 
alised, albeit of far higher resolution than previously doc¬ 
umented cosmological-scale simulations whose datasets en¬ 
compass this scenario, motivating the following discussion 
with proposals for future refinements to our model. 


4.1 A comparison of group tidal mechanisms with 
previous studies 

Our study is novel with respect to the adopted satel¬ 
lite/host masses, but is otherwise consistent with previ¬ 
ous numerical studies that demonstrate how satellites on 
eccentric orbits of a massive neighbour with weak dy¬ 
namical friction experience a succession of impulsive tidal 
shocks (due to the group tides) which heat the stel¬ 
lar disc (Gnedin, Hernquist & Ostriker 1999). This mech¬ 
anism will depend on the orbital pericentre radius/period 


(McGarthy et al. 2008; Villalobos et al. 2012), with those 
satellites that generally avoid the densest (central) regions 
of their respective environment, such as those modelled in 
this study, losing their disc (and rotational-support) least 
efficiently (Mastropietro et al. 2005). 

An integral part of this tidal transformation is the 
triggering of a bar instability (Kazantzidis et al. 2011). 
This is quite resolution-dependent, in light of the stochas¬ 
tic nature of intrinsic disc evolution in numerical models 
(Sellwood V Debattista 2009), so the robustness of our Iso¬ 
lated simulations to this instability over 10 Gyr provides 
some confidence in our modelling. We predict therefore the 
tidal triggering of a bar in a 10^ M© satellite on its most 
common infall trajectory (Villalobos et al. 2012). The ef¬ 
ficiency with which stellar bars are destroyed by the sub¬ 
sequent growth of a central mass concentration/angular 
momentum exchange from infalling gas is also resolution- 
dependent, but our results (a bar lifetime of order 1-2 Gyr) 
are consistent with higher-resolution studies which consider 
this problem in detail (Bournard, Gombes V Semelin 2005). 

Assuming therefore the slow strangulation of a 10^ 
M© group satellite, tidal heating facilitated by our shorter- 
lived bar will be therefore less efficient than that demon¬ 
strated in pure collisionless simulations, such as those of 
Kazantzidis et al. (2011) who assumed the early and com¬ 
plete removal of gas in their models of dSph formation from 
10^ M© progenitor discs. 

Our approach to modelling a range of harassment sce¬ 
narios by other satellites (Section 3.3) does not highlight any 
signihcant systematic variation to this evolution. A similar 
result was found by Villalobos, De Lucia V Murante (2014) 
for more massive satellites, although their explicit modelling 
of a live group halo introduced the concept that other satel¬ 
lites can indirectly mitigate satellite transformation by in 
fact modifying the mutual group potential and the associ¬ 
ated satellite orbits. 

We assert therefore that for a given infall trajectory, 
dwarf satellite transformation is to hrst order dependent on 
time resident within the group halo. The concurrent sup¬ 
pression of SF faciliated by group tidal mechanisms is ineffi¬ 
cient, but an inverse correlation in the dwarf quenched frac¬ 
tion with mass (Wheeler et al. 2014) is consistent with the 
earlier infah times for lower mass satellites (De Lucia et al. 
2012; Wetzel et al. 2013), which would experience a smaller 
dynamical friction timescale in the lower mass hosts to which 
they are preferentially accreted at higher z (e.g. Jiang et al. 
2008). 

This assumes however that we have properly charac¬ 
terised the tidal interaction history of a 10^ M© satel¬ 
lite i.e. the incidence of major mergers, which we do not 
model in detail but which can signihcantly disrupt the stel¬ 
lar system (and form a slow rotating system; Fig. 9) is 
low. If adopting the dehnition of a major merger utilised 
by Rodriguez-Gomez et al. (2015) (e.g. with mass ratio > 
0.25), and integrating their estimated merger rate for this 
stellar mass which they deduce from the tracking of sub¬ 
halo merger trees in the Ihustris cosmological-scale simula¬ 
tion, we hnd up to 25 percent of such systems, a minority, 
will be involved in a major merger since z = 2. This is 
consistent with most slow rotating systems, whose forma¬ 
tion is linked with major mergers, being more massive sys¬ 
tems than considered here 10io-5 Emsellem et al. 
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2011), and whose mergers rate are higher for all z. More¬ 
over, bright dwarves corresponding to our 10^ Mq sys¬ 
tem are least likely among all dwarves to have been pre- 
processed and accreted as a satellite to a group-mass host 
(e.g. Wetzel, Deason & Garrison-Kimmel 2014), a process 
during which mergers/strong tidal interactions are feasible 
(Gnedin 2003). 

The tidal model of the group presented in this study 
is similar to that of Bekki & Gouch (2011), who demon¬ 
strated the transformation of spirals to SOs in response to 
recent findings that SOs are preferentially formed in (cr < 750 
kms“^) groups (Willman et al. 2009; Just et al. 2010). This 
supports the role of tidal interactions in their formation, 
as opposed to the hydrodynamical ISM-IGM interactions 
which are expected to be dominant in higher-cr environ¬ 
ments. The long transition timescale conveyed by our satel¬ 
lites is consistent with the continuous growth in SO numbers 
since intermediate redshifts (Just et al. 2010), while the sig¬ 
nificant reduction in Ar and s (Fig. 9) is consistent with 
tidal mechanisms being essential in their formation (as op¬ 
posed to being simply faded spirals, e.g. van der Wei et al. 
2010; Maltby et al. 2015). 

We stress however that we do not reproduce the classic 
morphology of SO in our studies, one discrepancy being the 
lasting spiral/tidal arm structures in our Tidal simulations 
(in addition to the ongoing SF which we discuss further in 
Section 4.2). The inability of our group environment to con¬ 
sistently form SOs is qualitatively in agreement with their 
being typically more massive than the stellar mass range 
considered here, as implied by their luminosity functions 
(e.g. Kelvin et al. 2014). Furthermore, the SO number frac¬ 
tion peaks in group systems with velocity dispersions of 600- 
700 kms“^ (Just et al. 2010), which are nominally far more 
massive, and quite likely to have accreted from, poor 10^^ 
Mq groups with a ~ 150 kms“^ (McGee et al. 2009). 

Instead, our Tidal simulations yield systems which are, 
morphologically and kinematically, characteristic of dwarf 
SOs. Unlike SOs, these systems are distinguished from the 
dE classification in which they are often subsumed by a two- 
component fit to their light profile and, saliently, evidence of 
a disc including spiral structure (Sandage & Binggeli 1984; 
Aguerri et al. 2005). 

Our simulations also convey how tidal influences within 
a 10^^ Mq group will be accompanied with ongoing star- 
bursts that are increasingly concentrated in the disc centre 
due to gas infall facilitated by tidal torques, bar dynam¬ 
ics and an enhanced Q in the outer disc. Optically, this 
would replicate the blue cores found among 70 percent of 
the dE/dSO galaxies of the Ursa Major cluster by Pak et al. 
(2014), with the dSO subsample exhibiting the most blue g—i 
colours for a gven magnitude (in a range —18 < Mr < —9). 
Since this cluster possesses no significant evidence of ongo¬ 
ing ram pressure stripping/H I-deficiency, they hypothesize 
that these galaxies represent satellites in transition, whose 
evolution is primarily driven by tidal interactions. This is 
supported by a correlation in the frequency of dEs/dSOs 
with the dynamical state/virialisation of their respective 
subgroups (see also Tully & Trentham 2008), and a system¬ 
atically more blue colour compared with their counterparts 
in the Virgo cluster (the latter exhibiting strong evidence 
of ram pressure stripping; Ghung et al. 2007). Our earlier 
assertion that the bright dwarves modelled here are least 


efficiently transformed by group tidal mechanisms is also in- 
triguingly supported by a statistically-significant (at the 3a 
level) absence of dEs in a range —18 < Mr < —17 that 
closely corresponds to our considered stellar mass. 

4.2 A comparison of group quenching mechanisms 
with previous studies 

Our simulations indicate that the tidal mechanism acting 
within an idealised group (e.g. harassment, stripping) cannot 
consistently quench a bright dwarf satellite, in which case 
the environmental processes that dominate the formation of 
quiescent systems (Wetzel et al. 2013), in particular those in 
the dwarf mass regime (Geha et al. 2012), hinge upon the 
interaction of the satellite with the intra-group medium. A 
need to characterise its efficiency is highlighted by increasing 
evidence that 10^^“^^ M© hosts are the building blocks of 
larger environments, and where quenching mechanisms ap¬ 
pear to commence (McGee et al. 2009; Rafieferantsoa et al. 
2015). 

Our assumption of the instantaneous removal of a 10^ 
M© satellites’ hot gaseous corona by the group IGM is 
broadly consistent, albeit conservative, with respect to pre¬ 
vious studies which explicitly model this phenomenon and 
find it almost complete within a Gyr (Kawata V Mulchaey 
2008; Bekki 2009). This is supported observationally by 
the apparent suppression in satellites on their first infall to 
groups (Rasmussen, Ponman & Mulchaey 2006), and even 
up to several rvir from group centres (Rasmussen et al. 
2012). The suppression of SE by the shutting off of cold 
mode gas accretion predicted to contribute ^ 50 percent 
to the SER of bright dwarves (Sanchez Almeida et al. 2015) 
is supported further by recent cosmological simulations that 
highlight how the SE of galaxies is determined by this exter¬ 
nal gas source prior to infall (Rafieferantsoa et al. 2015). hol¬ 
lowing accretion to halos < 10^^ M©, Rafieferantsoa et al. 
demonstrate how the H I-richness of galaxies traces the 
sSER, implying the starvation of the satellite on a timescale 
of several Gyr. 

Various studies espouse a similar timescale during which 
the SER declines in an exponential manner since infall, 
with e-folding timescales inferred to be of order 1-3 Gyr 
(Wang et al. 2007; McGee et al. 2009; De Lucia et al. 2012) 
but which in practice can vary widely depending on respec¬ 
tive satellite and host masses (Bekki 2014). This lies broadly 
in accordance with the predominantly secular consumption 
of extant disc gas in our Tidal simulations irrespective of 
triggered starbursts (Eig. 6), which proceeds according to a 
SE model that matches a dependence on gas surface density 
for analogous metal-poor dwarf galaxies (Bigiel et al. 2008). 
The low SE efficiency with respect to e.g. the Kennicutt- 
Schmidt relation (Eig. 4) is thus demonstrated to match the 
slow transition to an inactive status.. 

An alternative two-stage mode of suppresion was re¬ 
cently introduced by Wetzel et al. (2013), who compare the 
SER of satellites with a statistical parametrization of central 
galaxy SERs in the SDSS DR 7 (with a completeness limit 
lying slightly above our stellar mass at 10^'^ M©). Their 
analysis indicates that the suppresion of SE in dwarf satel¬ 
lites can be delayed by more than 4 Gyr after infall (during 
which the build up of new stellar mass is substantial, of or¬ 
der 50 percent), followed by a rapid decline in SER with 
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timescale < 1 Gyr. Their conclusion that these satellites are 
accreted with enough cold gas to sustain the SFR prior to 
infall is not supported by our simulations, although we have 
adopted an initial gas fraction (fgas) implied by observations 
for this stellar mass (~ 1; Popping, Behroozi & Pee 2014). 
This may underestimate the total associated gas, which has 
recently demonstrated for infalling satellites of the Illustris 
simulation to be in the range fgas = 2 — 8 (Sales et al. 2015). 

We do note however a quantitative agreement in this 
two-stage satellite quenching with our ram pressure simula¬ 
tions (Section 3.4). Fig. 13 conveys a high sustained SFR (for 
several Gyr) in the Ram/z = 1/Overdensity simulation, be¬ 
fore the available gas is exhausted and/or the local ram pres¬ 
sure at pericentre passage is sufficient to effectively quench 
the satellite within 5 Gyr of first virial crossing. This trigger¬ 
ing of SF at the IGM-ISM interface, with simultaneous strip¬ 
ping in the outer disc, contributes to the concentration of SF 
at its centre where the restoring force is largest (Bekki 2014). 
A similar enhancement in SFR and a truncation in the 
Ha scalelength, compared to field counterparts, was noted 
in analyses of Virgo cluster spirals (Koopmann V Kenney 
2004; Koopmann, Haynes V Gatinella 2006). This phe¬ 
nomenon is not limited to dense environments, as re¬ 
vealed by the observation of both a starburst on the lead¬ 
ing edge and rapid gas loss via ram pressure stripping in 
the trailing tail of a late-type satellite of a poor group 
by (Rasmussen, Ponman V Mulchaey 2006). On the other 
hand, an analysis by Rasmussen et al. (2012) of SF rates 
among the constituents of 23 local galaxy groups finds only 
1 to 10 percent of group galaxies exhibiting evidence of 
starbursts; instead, the net effect on SF, with respect to 
their field sample, is a suppression by ~40 percent at group¬ 
centric distances up to 2rvir. 

A sustained triggering of SF is limited to a model 
in which we assume a locally overdense IGM, which re¬ 
cent large-scale hydrodynamical simulations predict infalling 
satellites to preferentially lie within, with the corresponding 
quenching timescale in group hosts less than the first orbital 
period (Bahe & McGarthy 2015; Gen, Pop V Bahcall 2014). 
However, their results, implying the complete quenching of 
the population of satellites of M© groups 

within 4 Gyr contrast sharply with those of Wheeler et al. 
(2014), who adopt mock observations from the Millenium- 
H simulation and establish the long quenching timescales 
(^6 — 7 Gyr) that motivated this paper. At present, we em¬ 
phasise that these cosmological simulations utilise baryonic 
particle masses of order 10^h“^ M© (corresponding to a res¬ 
olution of only ~100 particles for a stellar mass of 10® M©). 
Nonetheless, we identify some promising agreement with our 
better resolved but idealised group infall model that encour¬ 
ages further refinement. 
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